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ABSTRACT 

Using annual time series data on neonatal deaths in Zimbabwe from 1966 to 2018, we model and forecast 
number of neonatal deaths over the next 25 years using the Box — Jenkins ARIMA technique. Diagnostic tests such 
as the ADF tests show that Neonatal Deaths (ND) series is I (2). Based on the AIC, the study presents the ARIMA 
(8, 2, 0) model as the optimal model. The diagnostic tests further indicate that the presented model is stable and its 
residuals are stationary in levels. The results of the study reveal that the numbers of neonatal deaths per year are 
expected to decline sharply in the next 25 years. In order to keep on reducing neonatal deaths in Zimbabwe, the 
study offered a four-fold policy prescription. 
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1. INTRODUCTION 
Neonatal death can be defined as the number of neonates dying before reaching 28 days of age (Usman et al. 

2019). The first 2 days after birth account for over 50% neonatal deaths, while the first week of life accounts for 
over 75% of all neonatal deaths (Carlo & Travers, 2016). Actually, the risk of neonatal death is highest in the first 
24 hours of life (Nouri et al., 2013). In fact, 2.6 million children died in the first month of life in 2016 — nearly 
7000 newborn deaths every day — most of which occurred in the first week, with about | million dying on the first 
day and close to 1 million dying within the next 6 days (UNICEF, 2017). The major causes of neonatal deaths are 
birth asphyxia, prematurity, sepsis as well as congenital malformation (Carlo & Travers, 2016). Thus, neonatal 
deaths are an indicator of healthcare systems in a country (Babaei et al. 2018) in the sense that neonatal deaths 
reveal the health of children and development of the economy and culture of a country or region (Chengye, 2012). 
Interestingly, neonatal deaths are preventable (Tachiwenyika et al., 2011). In order to enhance the prevention of 
neonatal deaths, modeling and forecasting neonatal deaths is critical, especially in developing countries such as 
Zimbabwe where neonatal deaths are still prevalent in large numbers. Therefore, this paper, will go a long way in 
uncovering the dynamics of neonatal deaths in Zimbabwe and consequently shed more light on health policy 
formulation. 
1.22 OBJECTIVES OF THE STUDY 

1. To investigate the years during which neonatal deaths peaked in Zimbabwe. 

il. To forecast the number of neonatal deaths for the out-of sample period. 

i. To examine the pattern of neonatal deaths for the out-of-sample period. 


1.3 RELEVANCE OF THE STUDY 

Neonatal death is still a major public health problem worldwide and accounts for more than 60% of 
newborn deaths before their first birthday (UNICEF, 2008). Of the world’s 7.7 million deaths in those aged 
younger than 5 years, 3.1 million are neonatal deaths (Rajaratnam et al. 2010). Approximately 99% of these 
neonatal deaths occur in low and middle — income countries, mostly in sub-Saharan Africa (Lawn et al. 2005) 
including Zimbabwe which continues to bear a heavy burden of neonatal mortality (Ministry of Health and Child 
Care, 2007). This study seeks to examine and forecast the number of neonatal deaths in Zimbabwe. In order to 
reduce the numbers of neonatal deaths to zero, there is need for reliable forecasts that will act as a guiding tool 
for policy makers in the health sector; hence, the need for this study. 


2. LITERATURE REVIEW 
Sarpong (2013) modeled and forecasted maternal mortality ratio (MMR) at the Okomfo Anokye 
Teaching Hospital in Kumasi, Ghana, from the year 2000 to 2010; using ARIMA models and found out that the 
ARIMA (1, 0, 2) model was optimal for forecasting quarterly MMR at Okomfo Anokye Teaching Hospital. Ezeh 
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et al. (2014) analyzed the determinants of neonatal mortality in Nigeria using the Cox Regression model and found 
out that a higher birth order of newborns with a short birth interval of less or equal to 2 years and newborns with 
a higher birth order with a longer birth interval of greater than 2 years were significantly associated with neonatal 
mortality. 

Nyoni (2019) modeled and forecasted maternal deaths in Zimbabwe using annual time series data 
covering the period 1990 — 2015 and applied the Box-Jenkins ARIMA models and basically found out that in the 
next decade (2016-2025), maternal deaths will increase. In another Zimbabwean study, Chaibva et al. (2019) 
analyzed stillbirths and neonatal deaths in Mutare district: the study conducted a retrospective review of 346 
patient records, of women who delivered at Sakubva Hospital and those reffered for Mutare district facilities to 
Mutare Provincial Hospital, between January and June 2014 and then used descriptive statistics to explore the 
contributors to stillbirths and neonatal deaths in Mutare. Their results basically show that of the 346 women, 
15.6% (i.e. 54) experienced an adverse pregnancy outcome (stillbirth or neonatal death). Their results also indicate 
that contributing factors to adverse pregnancy outcomes included birthweight, gestational age, delivery 
complications and delivery methods. In yet another, most recent Zimbabwean study, Nyoni & Nyoni (2020) 
analyzed monthly time series data on neonatal death cases at Chitungwiza Central Hospital (CCH) from January 
2013 to December 2018 using Box-Jenkins SARIMA models and found out that there will be a slow but steady 
decrease in neonatal deaths at CCH over the out-of-sample period, that is, January 2019 to December 2020. 

Mishra et al. (2019) forecasted Infant Mortality Rates (IMR) in India using ARIMA models. The forecast 
of the sample period (1971-2016) indicated accuracy by the selected ARIMA (2, 1, 1) model. The post sample 
forecast with the ARIMA (2, 1, 1) model revealed a decreasing trend of IMR (2017-2025). The forecast IMR for 
2025 was found to be 15/1000 live births. Khan et al. (2019) modeled and forecasted IMR of Asian countries 
using the log-log regression and ARIMA models and found out that there was a negative correlation between IMR 
and GDP (PPP). Secondary data of IMR and GDP (PPP) from 1980 to 2015 was analyzed and forecast was done 
from 2016 to 2025: the AR (1) model was found for all countries except Japan and Nepal for which the ARIMA 
(1, 1, 1) model was found suitable. Usman et al. (2019) analyzed the incidence of the rate of neonatal mortality in 
Nigeria using ARIMA models. Their trend plot of the incidence indicated that there was a steady decrease in the 
incidence rate over the years. The ARIMA (1, 1, 1) model was found to be the optimal model. The time series 
analysis also revealed that the neonatal mortality rate has reduced by 17.8% from 51.7% in the year 1990 to 33.9% 
in the year 2017. This paper follows the leads of Usman et al. (2019) and is the first country-specific study which 
has forecasted neonatal deaths in Zimbabwe. 


3. MATERIALS & METHODS 
ARIMA Models 
Autoregressive Integrated Moving Average (ARIMA) models deliver more accurate forecasts than 
econometric techniques (Song et al., 2003b). In fact, ARIMA models perform better than multivariate models in 
forecasting (du Preez & Witt, 2003). ARIMA models were developed by Box & Jenkins (1970) and their approach 
of identification, estimation and diagnostics 1s based on the principle of parsimony (Asteriou & Hall, 2007). The 
generalized ARIMA (p, d, q) model can be represented by a backward shift operator as: 


O(B)G. = B)°ND p= OCB) fle ee diese ssa ativeirsuccievss ued aderate dune var vente dan athasieue iver steams [1] 


Where the autoregressive (AR) and moving average (MA) characteristic operators are: 


CL Be BP tdci pee eos eb tae teeters eee ei | 
OB) CG Be Be eae eee ane eae) 
and 
(1 — B)INDy, = AGNDg ose cee cee cee cee cee cee ee cee cee cee ce cee see aes cee cessor ees cassie at ens sie cian ane ceee [4] 


Where @the parameter estimate of the autoregressive component is, 0 is the parameter estimate of the 
moving average component, A 1s the difference operator, d is the difference, B is the backshift operator and [; is 
the disturbance term. 


The Box — Jenkins Methodology 

The first step towards model selection is to difference the series in order to achieve stationarity. Once 
this process is over, the researcher will then examine the correlogram in order to decide on the appropriate orders 
of the AR and MA components. It is important to highlight the fact that this procedure (of choosing the AR and 
MA components) is biased towards the use of personal judgement because there are no clear — cut rules on how 
to decide on the appropriate AR and MA components. Therefore, experience plays a pivotal role in this regard. 
The next step is the estimation of the tentative model, after which diagnostic testing shall follow. Diagnostic 
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checking is usually done by generating the set of residuals and testing whether they satisfy the characteristics of 
a white noise process. If not, there would be need for model re — specification and repetition of the same process; 
this time from the second stage. The process may go on and on until an appropriate model is identified (Nyoni, 
2018c). 
Data Collection 

This study is based on 53 observations of annual total Neonatal Deaths (ND) in Zimbabwe. All the data 
was gathered from the World Bank online database. 


Diagnostic Tests & Model Evaluation 
Stationarity Tests: Graphical Analysis 





Figure 1 


Figure | above indicates that the ND series is not stationary since it follows a particular trend over the 
period 1966 to 2018. This basically implies that the mean and varience of the ND series is changing over time. 
Between 1966 and 1980, neonatal deaths were on the rise in Zimbabwe (then Rhodesia). This could be attributed 
to the liberation war (between black majority and white minority) that was taking place in Rhodesia. Soon after 
Zimbabwe’s independence, the country inherited a health system which was well functioning and given there was 
political stability; neonatal deaths dropped significantly from as high as 10869 deaths in 1980 to as low as 8455 
in 1993. The disastrous macroeconomic reforms over the period 1990 — 2000, largely contributed to poor 
performance of the health sector and hence neonatal healthcare service delivery was worse off. The following 
“lost decade’, that is; 2000 to 2010 was a completely lost decade, as noted by Kanyenze et al. (2017), as it was 
characterized by gross macroeconomic mismanagement, hyperinflation and excessive unemployment. This period 
was a huge blow to the health sector in Zimbabwe and this could be an explanation as to why neonatal deaths had 
to sky-rocket over the period 2000 to 2010. Thereafter, the numbers of neonatal deaths started going down 
gradually. This could be attributed to macroeconomic stability that was largely brought about by the introduction 
of the United States Dollar (USD) as the official currency, following the rejection of the Zimbabwean dollar which 
had lost value. When the economy is performing, the government and its partners are able to mobilize resources 
for the health sector and this improves health service delivery. When the economy is not performing, health 
workers migrate to greener pastures just like what happened during the lost decade. The government is usually 
not able to capacitate and renovate existing healthcare facilities if the economy is not performing. In order to 
determine the order of integration of the ND series shown above, the study will employ correlogram analyses 
along with the Augmented-Dickey-Fuller (ADF) test. 
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The Correlogram in Levels 


ACF for ND 


+- 1.06/7° 0.5 


PACF for ND 


+- 1.06/T7°0.5 





Figure 2 
The ADF Test 
Table 1: Levels-intercept 


po 2.923780 | 5% | Stationary 
pd 2599925 | @10% | Stationary 


Table 2: Levels-trend & intercept 


pO 3.506374 | 5% | Stationary 
po 3.183002 | @ 10% | Stationary 


Table 3: without intercept and trend & intercept 


ADE Statistic Probability Critical Values 
0.305604 0.7702 2.613010 











| 
po 1.947665 
po 1.612573 @10% 
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The Correlogram (at 1‘ Differences) 


ACF for d_ND 


+: 1594 2.5 —= 


PACF for d_ND 


+: 1964 0.5 —= 





Figure 3 


Table 4: 1‘ Difference-intercept 


3.571310 
po 2.922449 | @S% | Stationary 
po 2.599224 | @10%__| Stationary 


Table 5: 1‘ Difference-trend & intercept 


po 3.504330 | @5% | Notsstationary 
po 3181826 | @ 10% __| Not stationary 








Table 6: 1‘ Difference-without intercept and trend & intercept 


| Variable —_| ADF Statistic Probability Critical Values 
|ND__| -3.071985 0.0028 2.613010 
Ss 1.947665 





po 1.612573 @10% 


Figures above, that is; 2 and 3 and tables above, that is; 1 to 6 show that the ND series is not stationary in levels 
and even after taking first differences. 
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The Correlogram in (2" Differences) 


ACF for d_d_ND 


+- 1.06/10. 5 —= 


8 


lag 


PACF for d_d_ND 


+- 1.06/70. 5 —= 





Figure 4 


Table 7: 2™ Difference-intercept 


po 2.923780 | @5% (| Stationary 
po 2.599925 | @10% | Stationary 





Table 8: 2™ Difference-trend & intercept 


| Variable —_| ADF Statistic Probability Critical Values 
PND -3.551269 0.0452 4.161144 
po 3.506374 





3.183002 @10% 


Table 9: 2™ Difference-without intercept and trend & intercept 


ADF Statistic Probability Critical Values 
“3.656290 0.0005 2.614029 





po 1.947816 
pf 1.612492 @10% 


Figure 4 and tables 7 — 9 illustrate that the ND series is I (2). 
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Evaluation of ARIMA models (without a constant) 
Table 10: Evaluation of ARIMA Models 


mic fu [we [waz [RMSE [WARE 
28.506 


A model with a lower AIC value is better than the one with a higher AIC value (Nyoni, 2018b) Similarly, 
the U statistic can be used to find a better model in the sense that it must lie between O and 1, of which the closer 
it is to O, the better the forecast method (Nyoni, 2018a). In this paper, only the AIC is used to select the optimal 
model. Therefore, the ARIMA (8, 2, 0) model is chosen. 





Residual & Stability Tests 
ADF Tests of the Residuals of the ARIMA (8, 2, 0) Model 
Table 11: Levels-intercept 


Ro | 6.595173 0.0000 | 3.596616 | @1% | Stationary 
po 2.933158 | 5% | Stationary 
po 2.604867 | @ 10% _| Stationary 


Table 12: Levels-trend & intercept 


Ro | 6.521314 0.0000 | 4.192337 | @1% | Stationary 
po 3.520787 | @5%__| Stationary 
po 3191277 | @ 10% | Stationary 


Table 13: without intercept and trend & intercept 


ADE Statistic Probability Critical Values 
-R_ | -6.651100 0.0000 “2.621185 








-1.948886 Stationary 
1.61192 @ 10% 


Tables 11 — 13 indicate that the residuals of the chosen optimal model, the ARIMA (8, 2, 0) model; are stationary. 





Correlogram of the Residuals of the ARIMA (8, 2, 0) Model 
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Residual ACF 


+-)9106/T*0.5 


Residual PACF 


+- 91 06/T*0.5 





Figure 5: Correlogram of the Residuals 


Figure 5 indicates that the estimated model is adequate since ACF and PACE lags are quite short and 
within the bands. This implies that the no autocorrelation assumption is not violated in this study. 


Test for Normality of Residuals 


Test statistic for normality: uhat1 
Chi-square(2) = 0.205 [0.9024] N(-1.2423,27.985) 





Figure 6: Normality Test 
Since the p-value, that is; [0.9024] is statistically insignificant, it implies that the residuals are normally 
distributed, hence the validity of the normality assumption. 
Stability Test of the ARIMA (8, 2, 0) Model 
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Inverse Roots of AR/MA Polynomial(s) 





Figure 7: Inverse Roots 
Since all the AR roots lie inside the unit circle, it implies that the estimated ARIMA process is 
(covariance) stationary; thus confirming that the ARIMA (8, 2, 0) model is indeed stable and suitable for 
forecasting annual neonatal deaths in Zimbabwe. 


4. FINDINGS 
Descriptive Statistics 
Table 14: Descriptive Statistics 
Statistic 


Statistic 


As shown above, the mean is positive, 1.e. 9927. This means that the average number of neonatal deaths 
over the study period is 9927 deaths per annum. The minimum number of neonatal deaths over the study period 
is 7361 deaths and this was recorded in 1961 while the maximum number of neonatal deaths is 13169 deaths and 
this was recorded in 2010. The skewness is 0.51971 and the most important characteristic is that it is positive, 
meaning that the ND series is positively skewed and non-symmetric. Excess kurtosis is -0.48186; showing that 
the ND series is not normally distributed. 





Results Presentation 
Table 15: Main Results 
ARIMA (8, 2, 0) Model: 
A*ND, = 1.3784A7ND,_, — 0.410827A2ND,_, — 0.00721141A7ND,_3 — 4.29788A2ND,_, 
+ 0.556302A*ND,_; — 0.628168A7ND,_, + 0.648328A2ND,_- 
= 0.35054 20ND sii vie ai niareavaenrcwnenas .. [5] 


Standard Error 
1.37840 0.130935 10.53 0.0000" 


Be | 0.628168 | 0.232070 2.707 | 0.00608 


Forecast Graph 
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Figure 8: Forecast Graph — In & Out-of-Sample Forecasts 


Table 18: Tabulated Out-of-Sample Forecasts 


2025 
2026 
2027 
2028 
2029 
2030 
2031 
2032 
2033 
2034 
2035 
2036 
2087 
2058 
2039 


Predicted ND 
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Figure 9: Graphical Analysis of Out-of-Sample Forecasts 
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Table 15 shows the main results of the ARIMA (8, 2, 0) model. Figure 8 and 9 as well as table 18 are 
out-of-sample forecasts of the ARIMA (8, 2, 0) model. As clearly shown in figure 9, the number of neonatal 
deaths per year, over the out-of-sample period, show a sharply downwards trend. This is encouraging and 
commendable, for a developing country like Zimbabwe. These results are consistent with Nyoni & Nyoni (2020). 


Policy Implications 

1. The government of Zimbabwe should continue to intensify training programs in resuscitation and in essential 
newborn care in order to maintain low levels of and or eradicate neonatal deaths. 

ii. The government of Zimbabwe should work towards improving access to healthcare services through out the 
whole country. 

ii. The government of Zimbabwe should also work toward capacity building in virtually all public health 
institutions in the country to ensure that comprehensive neonatal care services are offered country-wide. 

iv. The government of Zimbabwe should encourage and promote consistent home visits by community health 
workers, for neonatal care. 


5. CONCLUSION 
The study shows that the ARIMA (8, 2, 0) model is not only stable but also the most suitable model to 
forecast neonatal deaths in Zimbabwe for the next 25 years. The model predicts a sharp decrease in neonatal 
deaths in Zimbabwe. Such a trend should be maintained and in this regard, a four-fold policy prescription has 
been offered. These findings are essential for the government of Zimbabwe, especially when it comes to long- 
term planning with regards to neonatal care in the country. 


REFERENCES 


1. Asteriou, D. & Hall, S. G. (2007). Applied Econometrics: a modern approach, Revised Edition, Palgrave 
MacMillan, New York. 

2. Babaei, H., Dehghan, M., & Pirkashani, L. M. (2018). Study of Causes of Neonatal Mortality and Its Related 
Factors in the Neonatal Intensive Care Unit of Iman Reza Hospital in Kermanshah during (2014 - 2016), 
International Journal of Pediatrics, 6 (5): 7641 — 7649. 

3. Box, G. E. P., & Jenkins, G. M. (1970). Time Series Analysis: Forecasting and Control, Holden Day, San 
Francisco. 

4. Carlo, W. A., & Travers, C. P. (2016). Maternal and Neonatal Mortality: Time To Act, Journal of Pediatrics, 
92 (6): 543 — 545. 


Volume 3, Issue VII, July 2020 | 49 


| ©) | Aad hd dEd e-ISSN : 2620 3502 
f ~~ International Journal on Integrate ucation p-ISSN : 2615 3785 


Pen 


5. Chaibva, B. V., Olorunju, S., Nyadundu, S., & Beke, A. (2019). Adverse Pregnancy Outcomes “Stillbirth and 
Early Neonatal Deaths” in Mutare District, Zimbabwe (2014): A Descriptive Study, BMC Pregnancy and 
Childbirth, 19 (86): 1 — 7. 

6. Chengye, J. (2012). Child and Adolescent Health, People’s Medical Publishing House, Beijing. 

7. du Preez, J. & Witt, S. F. (2003). Univariate and multivariate time series forecasting: An application to 
tourism demand, International Journal of Forecasting, 19: 435 — 451. 

8. Ezeh, O. K., Agho, K. E., Dibley, M. J., Hall, J., & Page, A. N. (2014). Determinants of Neonatal Mortality 
in Nigeria: Evidence From the 2008 Demographic and Health Survey, BMC Public Health, 14: 521 — 531. 

9. Kanyenze, G., Chitambara, P., & Tyson, J. (2017). The Outlook For The Zimbabwean Economy, Supporting 
Economic Transformation (SET), Harare. 

10. Khan, M. S., Fatima, S., Zia, S. S., Hussain, E., Faraz, T. R., & Khalid, F. (2019). Modeling and Forecasting 
Infant Mortality Rates of Asian Countries in the Perspective of GDP (PPP), International Journal of Scientific 
and Engineering Research, 10 (3): 18 — 23. 

11. Lawn, J. E., Cousens, S., & Zupan, J. (2005). Neonatal Survival 1:4 Million Deaths: When? Where? Where? 
Why” Neonatal Survival Series Paper 1, Lancet, 365: 891 — 900. 

12. Ministry of Health and Child Care (2007). The Zimbabwe National Maternal and Neonatal Health Road Map 
(2007-2015), Government of Zimbabwe, Harare. 

13. Mishra, A. K., Sahanaa, C., & Manikandan, M. (2019). Forecasting Indian Infant Mortality Rate: An 
Application of Autoregressive Integrated Moving Average Model, Journal of Family and Community 
Medicine, 26: 123 — 126. 

14. Nouri, A., Barati, L., Qhezelsofly, F., & Niazi, S. (2013). Causes of Infant Mortality in Kalaleh City During 
2004 — 2012, Hakim Jorjani Journal, 1 (2): 2 — 37. 

15. Nyoni, S. P., & Nyoni, T. (2020). ARIMA Modeling of Neonatal Mortality in Chitungwiza Central Hospital, 
International Journal of Multidisciplinary Research IJMR), 6 (2): 189 — 196. 

16. Nyoni, T (2018b). Modeling and Forecasting Inflation in Kenya: Recent Insights from ARIMA and GARCH 
analysis, Dimorian Review, 5 (6): 16 — 40. 

17. Nyon, T. (2018a). Modeling and Forecasting Naira/USD Exchange Rate in Nigeria: A Box-Jenkins ARIMA 
Approach, MPRA Paper No. 88622, University Library of Munich, Munich. 

18. Nyon, T. (2018c). Box — Jenkins ARIMA Approach to Predicting net FDI inflows in Zimbabwe, MPRA 
Paper No. 87737, University Library of Munich, Munich. 

19. Nyoni, T. (2019). Maternal Deaths in Zimbabwe: Is it a Crime to be a Woman in Zimbabwe? MPRA Paper 
No. 96789, University Library of Munich, Munich. 

20. Rajaratnam, J. K., Marcus, J. R., & Flaxman, A. D. (2010). Neonatal, postnatal, childhood and under-5 
mortality for 187 countries, 1970 — 2010: A Systematic Analysis of Progress Towards Millennium 
Development Goal 4, Lancet, 375: 1988 — 2008. 

21. Sarpong, S. A. (2013). Modeling and Forecasting Maternal Mortality; An Application of ARIMA Models, 
International Journal of Applied Science and Technology, 3 (1): 19 — 28. 

22. Song, H., Witt, S. F. & Jensen, T. C. (2003b). Tourism forecasting: accuracy of alternative econometric 
models, International Journal of Forecasting, 19: 123 — 141. 

23. Tachiwenyika, E., Gombe, N., Shambira, G., Chadambuka, A., Tshimanga, M., & Zizhou, S. (2011). 
Determinants of Perinatal Mortality in Marondera District, Mashonaland East Province of Zimbabwe, 2009: 
a Case Control Study, Pan African Medical Journal, pp: 1 — 8. 

24. UNICEF (2008). The State of the World’s Children, Child Survival-UNICEF. 

25. UNICEF (2017). The neonatal period is the most vulnerable time for a child. http://data.unicef.org/child- 
mortality/neonatal.html ; accessed 29/01/2020. 

26. Usman, A., Sulaiman, M. A., & Abubakar, I. (2019). Trend of Neonatal Mortality In Nigeria From 1990 to 
2017 Using Time Series Analysis, Journal of Applied Sciences and Environmental Management, 23 (5): 865 
— 869. 


Volume 3, Issue VII, July 2020 | 50 


